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Abstract. We report an unexpectedly strong and complex chemi¬ 
cal bonding ofrare-gas atoms to neutral gold clusters. The bonding 
features are consistently reproduced at different levels of approx¬ 
imation within density-functional theory and beyond: from GGA, 
through hybrid and double-hybrid functionals, up to renormalized 
second-order perturbation theory. The main finding is that the 
adsorption of Ar, Kr, and Xe reduces electron-electron repulsion 
within gold dimer, causing strengthening of the Au-Au bond. Dif¬ 
ferently from the dimer, the rare-gas adsorption effects on the gold 
trimeFs geometry and vibrational frequencies are mainly due to 
electron occupation of the trimer’s lowest unoccupied molecular 
orbital. For the trimer, the theoretical results are also consistent 
with far-infrared multiple photon dissociation experiments. 

Binding of rare-gas (RG) atoms to molecules and metal clus¬ 
ters has been studied intensively in the past years Due to their 
very stable closed-shell electronic configuration, when interacting 
with neutral species and clusters, RG atoms are expected to interact 
via dispersion forces and polarization by multipole moments of the 
clusters. 

However, a closer look to the binding of RG atoms to Au clusters 
reveals a much more complex nature of the interaction.^i^ In most 
of the previous studies, electrostatic effects are expected to play a 
major role in the interaction between the RG atom(s) and Au, since 
the Au atom is either charged or is a part of a polar molecule (with a 
large dipole moment). Reports on theoretical analysis of the bind¬ 
ing of RG atoms to bare neutral Aum clusters are very scarce.— 
We find that the present understanding of RG-Aum interaction is 
not conclusive due to the sensitivity of RG-Aum electronic struc¬ 
ture to the level of theory. Therefore, in this work we analyze 
bonding between RG atoms (Ne, Ar, Kr, and Xe) and the small¬ 
est Aum clusters (M = 2, 3) using a variety of different theoretical 
approaches. In fact, we believe that interaction of this type has been 
overlooked so far, at least for this class of system. Our analysis ex¬ 
plains several puzzling features of the observed spectra obtained 
with far-IR resonance-enhanced multiple photon dissociation (FIR- 
MPD) spectroscopy, which are reported in detail in Ref.— 

In Fig. [T]we show, in the upper half of each panel, the FIR-MPD 
spectra of AU 3 complexed with one or two Kr atoms, at r=100 K. 
The first striking aspect of the measured spectra is that the adsorp¬ 
tion of the second Kr changes the spectrum significantly. This fact 
suggests that the interaction of Kr with AU 3 cannot be treated as a 
perturbation. 

To explain these findings, we have calculated the geometry, elec¬ 
tronic structure, and vibrational frequencies of RG-Au complexes. 
The calculations were performed with FHI-aimsii program pack¬ 
age for an accurate all-electron description based on numeric atom- 
centered basis functions. The IR spectra at finite temperature were 
calculated by performing Born-Oppenheimer molecular dynamics 
simulations in the canonical ensemble and extracting from the tra¬ 
jectories the Fourier transform of the dipole-dipole autocorrelation 
function.- The computational details are given in— and SI. 

We first address the simplest system where we observe an un¬ 
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Figure 1. FIR-MPD^ (upper panels) and theoretical IR spectra— (lower 
panels) at T = 100 K. Theoretical harmonic spectra are also shown (as 
upside-down bars in the lower panels) for the two isomers of AusKr (top) 
and Au 3 Kr 2 (bottom). The yellow (dashed) lines refer to the obtuse-angled 
isomer, while the black lines refer to the acute-angled one. 


usually strong bond between Kr and a neutral, non-polar species: 
Kr adsorbed on Au 2 (note that adsorption of Kr on a single neutral 
gold atom leads to a weakly bound van-der-Waals, vdW, complex). 
Au 2 Kr’s geometry is linear and the Kr-Au equilibrium distance is 
found to be around 2.7 A, i.e., in the range of covalent bonding 
for gold atoms. Furthermore, upon adsorption of one Kr atom, the 
Au-Au distance is shortened by 0.3% and the IR spectrum becomes 
active (see Fig. I and Table I in SI), with a Au-Au stretching-mode 
frequency 5% blue-shifted with respect to the (IR-inactive) corre¬ 
sponding mode of the isolated dimer. We interpret these results 
as an indication of the strengthening of the Au-Au bond. To our 
surprise, both Mulliken and Hirshfeld analysis of the charge trans¬ 
fer (with both PBE and PBEO functionals) agree in assigning a net 
charge transfer from Kr to the dimer by 0.1 e~. Such transfer of 
negative charge from Kr to Au was also predicted in Ref.— on the 
basis of natural bond orbital analysis. This charge transfer is incon¬ 
sistent with the strengthening of the Au-Au bond; in fact, the Au- 
Au bond in the Au^ is weaker (i.e., longer and with softer stretch- 
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ing frequency, see SI) than in Au 2 . We find similar contradiction 
for Au 3 Kr^ complexes. Interestingly, the highest-frequency exper¬ 
imental peak (around 170 cm“^) in Fig. [T] blue-shifts upon second 
Kr adsorption. 

Geometries and binding energies obtained via PBE-i-vdW (i.e., 
the PBEJ^ functional, corrected for long-range vdW interactions 
via the Tkatchenko-Scheffier (TS) schemei^) are compared for Au 2 
and Au 2 Kr to PBEO-i-vdW and B3LYP hybrid functionals, M06- 
2X functional,— XYG3 doubly hybrid functional,— renormalized 
second-order perturbation theory (rPT2),— MP2, and CCSD(T). 
The molecular orbitals of AumKta^ clusters were also calculated 
with many-body perturbation theory methods GoWb@PBEO and 
self-consistent GW (sc-GW), as implemented in EHI-aims.i24iS 
CCSD(T) calculations were performed using Gaussian03 codeiS 
with the aug-cc-pVTZ-PP basis set^ and small-core relativistic 
pseudopotentials.— 

The calculated Au 2 Kr binding energies and geometries are sum¬ 
marized in Table[T] Compared to CCSD(T), we find that PBE+vdW 
performs reasonably well (3 % error in binding energies, and 1% in 
error bond lengths) and the XYG3 functional performs remarkably 
well. 

Table 1. Calculated properties of Au 2 and Au 2 Kr. The structure is relaxed 
at each level of theory, for the collinear geometry. The distance Jau-kt is the 
distance between Kr and the closest Au atom. The binding energy of Au 2 
is E^(Au 2 ) = E(Au 2 ) — 2E(Au). The adsorption energy of the Kr atom(s) 
onto the Au 2 is AE^(Au 2 Kr) = E(Au 2 Kr) — E(Au 2 ) — E(Kr), where E(.) is 
the total energy of the relaxed system. For all methods except PBE+vdW the 
energies are counterpoise corrected for the basis set superposition error. 
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PBE+vdW^ 

2.354 

2.509 

2.503 

2.728 

0.222 


PBEO+vdW^ 

2.050 

2.520 

2.514 

2.773 

0.161 


B3LYP 

2.055 

2.530 

2.528 

2.878 

0.070 


M06-2X 

1.432 

2.538 

2.528 

3.030 

0.130 


XYG3^ 

2.296 

2.486 

2.480 

2.740 

0.215 


rPT2@PBEa^ 

2.202 

2.500 

2.496 

2.785 

0.208 


MP2^ 

2.445 

2.429 

2.421 

2.620 

0.379 


CCSD(T)^ 

2.292 

2.484 

2.477 

2.685 

0.320 


Exp.i^ 

2.30 ±0.1 

2.470 

- 

- 

- 



All the employed methods agree in predicting shortening of the 
Au-Au bond upon Kr attachment. Consistently, the harmonic vi¬ 
brational frequency of the Au-Au stretch is increased (see Table I 
in SI). The binding energy of Kr to Au 2 is in the range 0.2-0.3 eV, 
which is much stronger than typical vdW binding (0.02 eV per Au- 
Kr pair, for this system, according to the TS scheme) 

The adsorption of other RG species was also considered. We 
find that Ne, Ar, Kr, and Xe bind by 0.02, 0.11, 0.22, 0.43 eV, 
respectively (with PBEO-i-vdW). In all cases the adsorption leads to 
a small but systematic shortening of the Au-Au bond and increase 
of the stretching frequency. 

To explain the above results, we analyze the changes in Kohn- 
Sham molecular orbitals upon Kr binding to Au 2 . As can be seen 
from Eig. [J] the orbitals of Kr atom are close in energy to the 
Oj 2 orbital of Au 2 , i.e., the linear combination of Au’s d^i orbitals. 
However, the overlap of these orbitals alone cannot lead to a cova¬ 
lent bonding, because both orbitals are fully occupied. This picture 
is not an artifact of a particular functional: in fact it is confirmed 
at the PBEO, GoWb@PBEO (Pigs. 4 and 5 in SI), and sc-GW (Pig. 
[3]) levels. Pollowing the RG series, we note (see Pig. 2 in SI) that 
the p valence orbital of Ar, Kr, and Xe, has an energy (with respect 
to the vacuum level) that is close to the energy of the Oj 2 orbital 
of the dimer. This is not the case for Ne and, as explained below, 
this is the reason for the much weaker (vdW only) bonding of Ne 
to Au 2 . 

Upon Kr adsorption, we find a reduction in the contribution of 
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Figure 2. Molecular-orbital picture based on atom- and angular- 
momentum-projected (PBE) electronic density of states, calculated for Kr 
bonding to Au 2 . The dashed lines link similar orbitals in different molecular 
arrangements. The labels 1 to 4 denote orbitals used to calculate the change 
in the partial electron density upon addition of two Kr to Au 2 , shown in Fig. 
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Figure 3. Evolution of self-consistent GW spectra from isolated Au 2 and 
Kr to relaxed Au 2 Kr. 


Au 5^ orbitals into the orbital of Au 2 (see Fig. [2]). By analyz¬ 
ing differences of orbital densities, we conclude that this reduction 
corresponds to a polarization of the ^ orbital towards Kr. The po¬ 
larization results in accumulation of the electronic charge between 
Kr and Au 2 . The polarisation does not occur for Ne because the 
mixing of Ne 2;? and Au2 g^^ is inhibited by the large difference 
in energies between these two orbitals. 

In order to quantify this accumulation, we integrated the total 
electron density difference between Au 2 Kr and Au 2 (at the geome¬ 
try of Au 2 Kr), within two planes passing through Au and Kr atoms 
and perpendicular to the molecular axis. With PBE, the integrated 
accumulation is 0.012^“. 

At the same time, the contribution of the cr|^ antibonding orbital 
to the G^^ antibonding orbital increases. This is particularly clear 
for the linear Au 2 Kr 2 complex, that has the same symmetry (Doo/?) 
as Au 2 . The adsorption energy of two Kr atoms, AE^, is 0.47 eV, 
slightly less than twice the adsorption energy of one Kr atom. The 
Au-Au bond is further strengthened: its length is decreased, and 
the Au-Au stretching frequency is increased (see Table I in SI). 
Fig.|4]shows that the increased mixing of the occupied and unoccu¬ 
pied antibonding orbitals of Au 2 results in the partial removal of the 
electrons and placing them to a more diffuse orbital. In other 
words, the electron-electron repulsion is reduced upon Kr attach¬ 
ment, allowing for more electron density to accumulate between 
the Au atoms. Overall, this redistribution of orbital electron densi¬ 
ties leads to an increase of the total density between the Au atoms. 
Indeed, using the same method as for the Au-Kr bond, we find an 
accumulation of 0.030^“ between Au atoms, consistent with the 
bond strengthening. We attribute the apparent charge transfer from 
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Kr to Au 2 , as given by both Mulliken and Hirshfeld analysis, to the 
diffuse nature of the Au 2 orbitals, resulting in a misleading assign¬ 
ment of charges to atoms. 

We now focus on AU 3 , that has two isomers,— both isosceles tri¬ 
angles, one obtuse-angled (oa), with a bond angle of about 140° and 
the other acute-angled (aa, bond angle of ~ 66 °). With PBE+vdW, 
the oa isomer is more stable by 0.12 eV. By adsorbing one or two 
Kr atoms to these two isomers, the two resulting structures have ap¬ 
proximately the same total energy (aa isomer less stable by 0.03 eV, 
with PBE-i-vdW). The final structures of Au 3 Kr and Au 3 Kr 2 are 
shown in the insets of Fig. [T] The comparison of theoretical and 
experimental spectra of AU 3 complexed with one and two Kr atoms 
has been already presented in Ref.— Here, we focus on the explana¬ 
tion of the blue-shift of the higher frequency peak, related to the aa 
isomer, when going from one to two adsorbed Kr (see Fig. [T] 

Attachment of one Kr to the aa isomer shortens the opposite Au- 
Au bond by 0.08 A, while adsorption of two Kr atoms shortens 
the adjacent Au-Au bond by 0.32 A. The corresponding adsorption 
energies are 0.23 and 0.34 eV. In order to explain the calculated 
binding energies, the changes in geometry, and the blue-shift of 
the highest-frequency mode upon Kr attachment to AU 3 , we per¬ 
form the analysis similar to the Au 2 Kr case. The highest-occupied 
molecular orbital (HOMO) of AU 3 is singly occupied. Its charac¬ 
ter is bonding for the two short bonds and antibonding for the long 
bond (opposite to the 66 ° angle). The lowest unoccupied molec¬ 
ular orbital (LUMO) of (relaxed) AU 3 is antibonding for the short 
bonds, and bonding for the other bond. In this case, the 4p orbital 
of single added Kr mixes with the LUMO of AU 3 , leading to the 
increased LUMO-J-manifold gap, decreased Au-Au distance for 
the long bond, and increased Au-Au distances for the equivalent 
short bonds (see Fig. 7 in SI). Kr 4/? orbital does not mix with 
HOMO of Au 3 by symmetry. Thus, contrary to Au 2 complexes, 
the bonding can be described as partial donation of a Kr electron 
pair to the AU 3 cluster. Addition of the second Kr atom leads to 
swapping of the LUMO and HOMO orbitals (so that the former 
LUMO becomes singly occupied), and a large reduction (—11%) 
of the bond distance of the longest bond in AU 3 (see Table I and 
Fig. 6 in the Suppl. Material). The vibrational analysis shows that 
it is the strengthening of this bond that leads to the blue shift of the 
observed high-frequency peak. 

In conclusion, we find that RG atoms bind covalently to small 
neutral Au clusters. For closed-shell clusters (Au 2 in particular), 
rather exotic features of this binding become prominent: it involves 
a redistribution of electrons among the d- and .^-derived antibond¬ 
ing orbitals, leading to the strengthening of the Au-Au bonds. In 
the open-shell complexes (such as Au 3 Kr), however, mainly the 
partial donation of electrons from Kr to the unoccupied orbitals of 
the Au cluster determines the bonding. Beside fully explaining fine 
details in the FIR-MPD of the gold trimer, our findings suggest that 
weakly adsorbing species can fine-tune the electronic properties of 
very small clusters, which may be of interest for the use of gold 


clusters as catalysts. 
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Figure 4. Contour plot of the linear combination of the density of 
the molecular orbitals (MO) of Au 2 Kr 2 and Au 2 , as labeled in Fig. 
12] (MO 3 +MO 4 )—(MO 1 +MO 2 ). The isolevels are equally spaced by 

0.002e-/A^ 
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